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Abstract 

A computational study of the ram accelerator, a 
ramjet-in-tube device for accelerating projectiles to ultra- 
high velocities, is presented. The analysis is carried out 
using a fully implicit TVD scheme that efficiently solves 
the Reynolds-averaged Navier-Stokes equations and the 
species continuity equations associated with a finite rate 
combustion model. Previous analyses of this concept have 
been based on inviscid assumptions. The present results 
indicate that viscous effects are of primary importance; 
in all the cases studied, shock-induced combustion always 
started in the boundary layer. The effects of Mach num- 
ber, mixture composition, pressure, and turbulence are 
investigated for various configurations. Two types of com- 
bustion processes, one stable and the other unstable, have 
been observed depending on the inflow conditions. In the 
unstable case, a detonation wave is formed, which propa- 
gates upstream and unstarts the ram accelerator. In the 
stable case, a solution that converges to steady-state is 
obtained, in which the combustion wave remains station- 
ary with respect to the ram accelerator projectile. The 
possibility of stabilizing the detonation wave by means 
of a backward facing step is also investigated. In addi- 
tion to these studies, two numerical techniques have been 
tested. These two techniques are vector extrapolation to 
accelerate convergence, and a diagonal formulation that 
eliminates the expense of inverting large block matrices 
that arise in chemically reacting flows. 

Introduction 

The ram accelerator is a ramjet-in-tube concept in 
which a shaped projectile can in principle be accelerated 
efficiently to velocities in excess of 10 km/s by means 
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of detonation waves or other shock-induced combustion 
modes. This concept, developed at the University of 
Washington 1,2 , can be scaled for projectile masses ranging 
from grams to hundreds of kilograms, and has the poten- 
tial for a number of applications, such as hypervelocity 
impact physics, direct launch to orbit of acceleration in- 
sensitive payloads, and hypersonic testing 3,4 . In the last 
application, the ram accelerator can be used as a hyper- 
sonic research facility for studying hypersonic aerodynam- 
ics and the supersonic combustion flowfields of interest to 
the National Aerospace Plane (NASP) program. An ex- 
perimental ram accelerator device, currently operating at 
the University of Washington 5,6 , has reached velocities in 
excess of 2.5 km/sec and Mach numbers as high as 8.4. 

Although several ram accelerator operation modes have 
been proposed 1 , the analysis in this paper will be centered 
on the “superdetonative” mode shown in Fig. 1. In or- 
der to operate in this combustion mode, the projectile 
must fly at superdetonative speeds, i.e., speeds above the 
Chapman-Jouguet detonation speed of the gas mixture. 
The gasdynamic principles that govern the flow and com- 
bustion processes in the superdetonative ram accelerator 
are similar to those related to hypersonic airbreathing 
propulsion systems, particularly to the oblique detonation 
wave engine 7 (ODWE). However, the device is operated 
in a different manner. 

In the superdetonative ram accelerator (Fig. 1), the 
centerbody is a projectile fired into a tube filled with 
a premixed gaseous fuel /oxidizer mixture. There is 
no propellant on board the projectile. Ignition of the 
fuel/oxidizer mixture is achieved by means of a series of 
shock waves that increase its temperature. When the 
ignition temperature is reached at a designed location, 
rapid chemical reactions release energy into the flowing 
stream. The energy addition will establish either a deto- 
nation wave or a shock-deflagration wave, depending pri- 
marily on the mixture composition, pressure and tube 
size. The combustion process creates a high pressure re- 
gion over the back of the projectile, producing a thrust 
force. The pressure, composition, chemical energy density 
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and speed of sound of the mixture can be controlled to op- 
timize the performance for a given flight condition. Since 
the fuel and oxidizer in the ram accelerator concept are 
premixed, the difficulties in obtaining rapid and complete 
mixing encountered by the ODWE (and the conventional 
scramjet) are circumvented. 

The performance of the superdetonative ram acceler- 
ator has been evaluated in the past by using one of the 
following two approaches: 1) simplified one-dimensional 
flow models, and 2) numerical simulations based on com- 
putational fluid dynamics (CFD) methods. 

Simplified one-dimensional models have been developed 
by Rom & Kivity 8 and Humphrey 9 . Their analysis is 
based on the jump conditions across the leading oblique 
shock and across the reflected (and assumed planar) det- 
onation wave. It is further assumed that the flow aft of 
the detonation is in a state of chemical equilibrium. Yip 
et. al. 10 included the effects of nonequilibrium chemistry 
by using a streamtube analysis close to the surface of the 
projectile. Their approach, however, decouples the chem- 
istry and the fluid dynamics and therefore does not take 
proper account of the effects of heat release on the shock 
wave structure. 

The main disadvantage of the above models is that they 
can be applied only to analyze the “nominal” operating 
condition. The ram accelerator is said to be operating at 
the nominal condition when the reflected shock wave from 
the tube wall intersects the shoulder (corner) of the pro- 
jectile. The expansion wave exactly cancels the reflected 
shock and no shock or expansion waves are transmitted 
downstream. For a given flight condition, optimum per- 
formance will be achieved by operating at the nominal 
design condition. Ram accelerator operation at the nom- 
inal condition, however, will occur only at the start of the 
projectile flight. As the projectile accelerates, the angle 
between the leading conical shock and the projectile will 
decrease, and the reflected shock (detonation) wave will 
impinge on the projectile at a location aft of the shoul- 
der, giving rise to a multiple shock-expansion pattern. 
In addition, reaching this nominal operating condition is 
difficult in practice, since at high Mach numbers the an- 
gle between the leading shock and the projectile is very 
small (about 3° at M = 8), and therefore the clearance 
between the projectile and the tube wall will have to be 
very small. The calculations of Humphrey 9 show that at 
a speed of 4 kin/sec, for example, the projectile diameter 
will have to be 92.4% of the tube diameter (compared to 
76% for a typical experimental setup) in order to attain 
the nominal operation condition. For the tube diameter 
of 38 mm considered, this implies that the clearance be- 
tween projectile and tube will be less than 1.5 mm, which 
is smaller than the boundary layer thickness predicted in 
the present study at similar conditions. 

Analysis of the superdetonative ram accelerator con- 


cept based on inviscid computational fluid dynamic meth- 
ods have been conducted by Bogdanoff &: Brackett 11 ’ 12 
and Yungster, et. al 13,14 . The CFD studies of Bogdanoff 
and Brackett were based on first and second order Go- 
dunov type differencing procedures. A global Arrhenius 
rate equation was used to model the combustion process, 
with the Arrhenius constants determined from experi- 
mental ignition delay studies. The results presented by 
Yungster, et. al., were obtained using the point implicit, 
total variation diminishing (TVD) MacCormack scheme. 
The analysis used a 7-species, 8-step reaction mechanism 
for hydrogen-oxygen combustion. 

The results presented in all the above studies have con- 
firmed the potential of the ram accelerator to efficiently 
accelerate large masses (up to hundreds of kilograms) to 
velocities in excess of 10 km/sec. In these studies, how- 
ever, the effects of viscosity, heat conduction and mass 
diffusion have all been neglected. This has been a neces- 
sary first step in determining the performance capabilities 
of the ram accelerator. However, recent studies conducted 
by the author on shock-wave/boundary layer interactions 
in premixed combustible gases 15 , indicate that viscous ef- 
fects can be important. In the present paper, numerical 
simulations of the superdetonative ram accelerator are 
obtained using the full Reynolds-averaged Navier-Stokes 
equations and a finite-rate combustion model. The effects 
of Mach number, mixture composition, pressure, and tur- 
bulence are investigated for various configurations. The 
numerical formulation used is discussed below, preceding 
the presentation of the computational results. 

Numerical Formulation 

A complete description of the numerical scheme used 
in this study has been previously presented in Ref. 15; 
therefore, it will only be discussed in general terms here. 

The Reynolds-averaged Navier-Stokes equations cou- 
pled with chemical nonequilibrium processes are consid- 
ered. The equations are written in nondimensional vari- 
ables and in generalized coordinates. The thermodynamic 
transport properties, such as the specific heat, thermal 
conductivity and viscosity for each species are determined 
by fourth-order polynomials of temperature. The thermal 
conductivity and viscosity of the mixture are calculated 
using Wilke’s mixing rule. The binary mass diffusivity 
between any two species is obtained using the Chapman- 
Enskog theory in conjunction with the Lennard-Jones in- 
termolecular potential functions. 

In the present study, a 7-species, 8-step reaction mech- 
anism for hydrogen-oxygen combustion is adopted. This 
model is a reduced reaction mechanism obtained from 
more complete models by the exclusion of the reactions 
involving and HO 2 , (which could be important in 

low temperature ignition studies). A complete descrip- 
tion of the reduced model and a discussion of its accuracy 
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and range of application can be found in Refs. 13-15. A 
more complete 9-species, 18-step combustion model has 
been recently implemented for testing several variations 
of the basic iterative scheme used in the present work (to 
be discussed below). 

The turbulent model adopted in the present study is 
the Baldwin-Lomax algebraic eddy viscosity model 25 and 
assumes constant turbulent Prandtl and Schmidt num- 
bers ( Pr t = Sc t = 0.9). This model is chosen for its 
simplicity and computational efficiency. 

The interactions between turbulence and chemistry, 
which enter into the numerical formulation through the 
source term W{, represent a very difficult problem. To ac- 
count for such interaction effects would require a closure 
method such as the probability density function (PDF) 
approach or a direct numerical simulation (DNS). Since 
effective PDF closure methods are not yet available and 
DNS methods are currently applicable only to relatively 
simple flows, the interactions between turbulence and 
chemistry are not considered in the present study. 

The fully coupled Navier-Stokes equations and species 
continuity equations are solved using a new fully implicit 
finite difference CFD code. The code employs an itera- 
tive method that is based on the lower-upper symmetric 
successive overrelaxation (LU-SSOR) implicit factoriza- 
tion scheme, and a second order symmetric total varia- 
tion diminishing (TVD) differencing scheme. To acceler- 
ate the convergence of the basic iterative procedure, this 
code can be combined with vector extrapolation methods, 
such as the Minimal Polynomial (MPE) and the Reduced 
Rank (RRE) Extrapolation. The extrapolation proce- 
dure solves a linear least squares problem and produces 
a sequence of approximations that, in general, has better 
convergence properties than the sequence obtained from 
the iterative scheme alone. A detailed description of these 
extrapolation techniques can be found in Refs. 16 and 17. 

Two different formulations of the LU-SSOR factoriza- 
tion scheme are currently implemented. In one formu- 
lation, the implicit operator includes the full Jacobian 
matrix of the chemical source term, leading to a precon- 
ditioner matrix of size n 3 xn,, where n, is the number 
of species; this matrix has to be inverted at every grid 
point. If the number of species considered is large, in- 
verting this preconditioner can be very expensive. There- 
fore, a second formulation has been introduced. In this 
technique, the Jacobian matrix is replaced by a diagonal 
matrix that is designed to approximate the time scaling 
effects obtained by using the full Jacobian. No matrix 
inversions are required in this formulation. Several di- 
agonalization methods were tested. The diagonalization 
method of Imlay, et. al 18 gave the best results. However, 
this method was found to be less efficient than the full 
Jacobian formulation. The details of the diagonalization 
method can be found in Refs. 18 and 19. 


Computational test of the numerical techniques 

Figure 2 shows the density residual history obtained 
with the two formulations (full Jacobian and Imlay’s di- 
agonalization) with and without extrapolation, for the 
case of a supersonic flow of a stoichiometric H 2 - air mix- 
ture past a compression corner. The chemical nonequi- 
librium processes are simulated by using a 9-species, 18- 
step finite-rate combustion model. When extrapolation 
is used, it is started after No iterations, and is imple- 
mented in the so called “cycling” mode, using a sequence 
of K max vectors obtained from the iterative scheme. The 
overhead in CPU time due to the use of extrapolation 
is very small (less than 1%) in the present case. The 
results indicate that savings of up to 40% in the overall 
computational work required to reach convergence can be 
realized by using RRE in combination with the basic it- 
erative scheme (using both formulations). Similar results 
are obtained with MPE. The results also indicate that 
for the present chemistry model, the diagonal formula- 
tion is less efficient, requiring approximately 45% more 
CPU time than the Full Jacobian formulation to reach 
convergence. 

Results 

When analizing shock-induced combustion phenomena, 
it is sometimes necessary to distinguish between a deto- 
nation wave and a shock-deflagration wave. In the lit- 
erature, different definitions have been used in the past 
for different problems. In the present study, the following 
terminology will be used: 

1) A shock-deflagration wave is defined as a shock-induced 
combustion process in which the interaction between 
shock and combustion is weak; the chemical reactions 
do not affect significantly the shock. In this process, the 
shock and the combustion front are essentially decoupled. 

2) A detonation wave is defined, following Pratt 20 , as a 
shock-induced combustion process in which “the combus- 
tion front follows so closely on the igniting shock wave 
that the two waves are fully pressure-coupled”. 

Strictly speaking, the second definition applies to a 
coupled shock-deflagration wave, since a “true” self- 
sustaining detonation wave has, in addition to being 
fully pressure-coupled, a characteristic nonsteady three- 
dimensional cell structure. Also, the first definition 
is sometimes applied to a decoupled shock-deflagration 
wave. 

In the present study, no distinction will be made be- 
tween a detonation wave (as defined above) and a self- 
sustaining detonation wave. Computations carried out 
on blunt projectiles fired into detonable gases have shown 
that the present numerical formulation can successfully 
reproduce the various shock-induced combustion regimes 
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observed in these flows. These regimes include decou- 
pled and coupled shock-deflagration waves, and combina- 
tions of self-sustaining oblique and overdriven detonation 
waves 13 . 

Turbulent Flow Calculations 

The computational studies of the superdetonative ram 
accelerator concept are carried out for configurations hav- 
ing dimensions similar to those of the experimental device 
presently operating at the University of Washington 1 ’ 5,6 . 
The first set of computations were conducted for the fol- 
lowing geometry and inflow conditions: 

Referring to Fig. 1, the projectile half angle cone is set 
to 9 = 14°, the tube diameter is d t = 3.0 cm, the maxi- 
mum projectile diameter is d v — 1.95 cm, and its length 
is L = 15 cm. The flow is assumed to be fully turbulent 
along the entire projectile. A constant projectile surface 
temperature T w — 600°K is also assumed. The gas mix- 
ture is stoichiometric H ^- air at a fill pressure of p^ — 1 
atm, and a temperature of Too = 300°K. 

Figure 3 shows the converged solutions obtained on a 
157x45 grid at three different Mach numbers. Nondimen- 
sional temperature contour lines (T/T a Q ) are plotted on 
the top half of the projectile, and nondimensional pres- 
sure color contours ( p/poo ) are plotted on the bottom 
half. For clarity, all contour plots of ram accelerator con- 
figurations are magnified in the vertical direction by a 
factor of 2. In Fig. 3a, the combustion process starts in 
the boundary layer at a location immediately behind the 
point where the reflected shock wave from the tube wall 
impinges on the projectile. The combustion that begins 
in the boundary layer propagates outwards and down- 
stream, and a shock-deflagration wave is established. 

For a higher Mach number flow ( M = 7.5), combustion 
begins prematurely in the boundary layer at the nose re- 
gion of the projectile, as shown in Fig. 3b. A very com- 
plex interaction between the shock-wave system and the 
chemically reacting boundary layer is observed. Com- 
plete combustion is achieved behind the shock being re- 
flected from the projectile surface. At a still higher Mach 
number, M = 8, combustion takes place along the en- 
tire boundary layer in the nose region of the projectile, 
as shown in Fig. 3c. Complete combustion is achieved in 
this case behind the first reflected shock wave from the 
tube wall. Note that in all of the three cases above, a high 
pressure region over the back of the projectile is created, 
and as a result, a positive thrust force is produced in all 
three cases. 

The pressure distribution along the projectile surface 
and tube wall for all three Mach numbers is shown in 
Fig. 4. Note again the higher pressures over the tail of 
the projectile as compared to the nose. For the M — 6.7 
and M — 7.5 cases the pressure increases in two stages, 


corresponding to the two shocks impinging on the projec- 
tile. Also note the small pressure disturbance created by 
the reacting boundary layer at the nose for the M — 7.5 
case. The maximum pressures are obtained on the tube 
wall for all three cases. The peak pressure observed for 
the M — 8 case has a value around forty times the fill 
pressure. 

Figure 5 shows the variation of y + at points nearest the 
wall, the skin friction coefficient, c/ = T/lp^U^, and the 
heat transfer coefficient , cj, = k /\ p x W 3 , . The value of 
the Reynolds number Re = pooUooL/ Poo for the M — 6.7, 
7.5 and 8.0 was 1.81 x 10 7 ,2.03 x 10 7 and 2.16 x 10 7 , 
respectively. Ideally, one should have y + — 0(1) every- 
where; this is the case for most of the constant diameter 
section of the ram accelerator. However, the boundary 
layer resolution is not good enough at the nose and tail 
sections resulting in too high y + values. Attempts to 
improve the resolution in these sections of the projectile 
resulted in computational cells with a very high aspect 
ratio (6x/6y > 1000) in the constant area part of the pro- 
jectile; this considerably slowed down or even prevented 
convergence. 

The skin friction coefficient plot, in Fig 5b, shows that 
combustion causes a reduction in the skin friction and cre- 
ates a high-temperature, low-density boundary layer that 
is more susceptible to separation than a similar nonreact- 
ing case. A separation bubble is observed for the M — 6.7 
case near the rear shoulder; a larger one occurs for the 
M = 8 case near the front shoulder. The heat transfer 
coefficient, shown in Fig. 5c, tends to follow the same 
qualitative variations as the pressure, decreasing across 
an expansion wave and increasing across a shock wave. 

Figure 6 shows the Mach number profile for the M = 
6.7 case at the x/L = 0.35 station which corresponds 
to the beginning of the constant diameter section of the 
projectile. This figure indicates that, for fully turbu- 
lent flow, the boundary layer thickness is approximately 
6/L « 0.012, or 6 « 1.8 mm. This fairly thick bound- 
ary layer suggest that operating the ram accelerator at 
nominal conditions may not be possible since the clear- 
ance between the projectile and the tube will have to be 
smaller than S. Further more, the boundary layer grows 
as the square of the Mach number, making the clearance 
problem even more severe as the projectile accelerates. 

A nondimensional pressure thrust, F, can be defined 
as: 


where F is the pressure thrust, is the fill pressure, and 
At is the tube area. The pressure thrust F is calculated 
by numerically integrating the pressure over the projectile 
surface. 

The nondimensional skin friction drag is similarly cal- 
culated by integrating the skin friction over the projectile 
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surface. The net thrust generated in the ram accelerator 
is calculated by subtracting the skin friction drag from 
the pressure thrust. Table 1 summarizes the results ob- 
tained with this configuration. It is important to point 
out that even for the cases where premature combustion 
occurs in the boundary layer at the frontal part of the 
projectile, a net positive thrust is still obtained. 

Also listed in Table 1 is the thrust pressure ratio, 4> t , 
defined as the net average drive pressure on the projec- 
tile (the net thrust divided by the maximum projectile 
cross-sectional area) divided by the maximum cycle pres- 
sure. This parameter provides a measure of the device’s 
launch capability versus the maximum pressure the pro- 
jectile and launch tube must withstand. 

The next calculation was conducted with a slightly dif- 
ferent configuration, having d t — 3.8 cm, and d p = 2.5 
cm, and using a mixture of 2 flj + O 2 4- 5 Ar. The ad- 
dition of a diluent, such as Argon, to the combustible 
mixture is typically done in order to change its speed of 
sound. By doing this, a different velocity range can be 
covered with the same projectile. 

In order to compare the results with those obtained 
in the previous case (at M = 6.7), the conditions were 
selected so that the maximum temperature in the bound- 
ary layer for nonreacting flow is about the same as in 
the previous M — 6.7 case (about 1200°K). Thus, the 
ignition process will occur at similar temperatures. The 
same peak temperature as the mixture was ob- 

tained for T w = 300°K, and M = 5.8. The Reynolds 
number was 1.72 x 10 7 in this case. For this mixture, the 
shock waves are stronger at a given Mach number than 
those obtained for the H^-a . ir mixture due to a higher 
value of 7 oo = ( c p /c v )oo (1.53 vs. 1.40). Starting from 
a nonreacting solution, the result of this calculation is 
shown in Fig 7 in the form of nondimensional temperature 
contours (top half) and Mach number contours (bottom 
half). The results are shown at various iteration stages 
starting from the nonreacting solution* . A completely 
different behavior is obtained in this case. The combus- 
tion process begins in the boundary layer (Fig. 7a), with 
no apparent separation. The combustion process seems 
to establish a detonation wave (Fig. 7b). Note the much 
steeper angle of this detonation wave compared to the 
shock-deflagration wave in Fig 3a. The detonation wave 
produces a separation of the boundary layer (Fig. 7c). 
The separation bubble continues to grow larger as the 

t The numerical method used in the present study is first or sec- 
ond order accurate in time for nonreacting flows. For chemically 
reacting flows, however, the equations become stiff and the nu- 
merical scheme solves the stiffness problem by preconditioning the 
equations, essentially rescaling them in time such that all chemi- 
cal and convective phenomena evolve on comparable pseudo-time 
scales. Therefore, the results shown in Fig. 7 are not strictly time 
accurate, but rather they show the evolution of the iteration process 
on a pseudo-time scale. 


detonation wave moves upstream (Figs. 7d and 7e), un- 
til what appears to be an overdriven detonation wave is 
formed (Fig. 7f). At this point, a large pressure will be 
produced over the frontal part of the projectile resulting 
in a net drag force. This phenomena is caused entirely 
by viscous interactions, since no such behavior was pre- 
viously observed in inviscid calculations 13,14 . 

Laminar Flow Calculations 

To determine the extent of the turbulence model effect 
on the flow and combustion processes in the ram accelera- 
tor, a series of calculations were performed assuming lam- 
inar flow. The geometry used in this case is d t = 3.6 cm 
and d p = 2.52 cm, and only the frontal part of the projec- 
tile was considered (0 < x/L < 0.61). Figure 8 shows the 
converged solution for a M — 6.5 flow in a stoichiometric 
H 2 -air mixture at — 1 atm, and T w = 600°K. The 
results show that a separation bubble is formed by the 
shock wave/boundary layer interaction, and that the com- 
bustion process remains confined to the boundary layer. 
Most of the combustion occurs downstream of the shock, 
but some combustion is seen to propagate upstream into 
the separation bubble. 

Next, poo is increased to 10 atm. The results for this 
case are shown in Fig. 9. Under these conditions, com- 
bustion upstream of the shock is enhanced by the higher 
pressure (Fig. 9b), and the temperature inside the sepa- 
ration bubble is much higher than for the poo — 1 atm 
case. The enhanced combustion enlarges the separation 
bubble and generates a secondary shock ahead of it. This 
shock wave is strong enough to ignite the main flow, and 
the shock/combustion wave moves upstream. 

The laminar flow results for a 2 H 2 + O 2 + 5 Ar mixture 
at poo — 1 atm and M — 6.5 are shown in Fig. 10. As 
for the turbulent case, the combustion process that starts 
in the boundary layer (Fig. 10a) evolves into a detona- 
tion wave (Fig. 10b). The combustion then propagates 
upstream through the boundary layer expanding the sep- 
aration bubble (Figs. 10c- lOe) and generating additional 
shocks as in the case of the if 2 -air mixture (Fig. 9). 

Thus, it seems that the unstable combustion mecha- 
nism is caused by a sequence of boundary layer sepa- 
ration, upstream propagation of the combustion process 
through the boundary layer, and generation of secondary 
shock/detonation waves. The process appears to be en- 
hanced by turbulence and by increasing the fill pressure. 

Laminar Flow Calculations with a Backward 
Facing Step 

A numerical study was conducted on a projectile having 
a backward facing step, with the purpose of testing if such 
a configuration could stabilize the combustion process, by 
serving as a barrier preventing the upstream propagation 
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of combustion. Fox this purpose, a two-block modifica- 
tion of the CFD code was developed and tested first on 
a simple parallel supersonic flow over a backward facing 
step. The calculations were compared with experimental 
results obtained by Donaldson 21 for a M = 3.5 laminar 
flow past a 1.9 cm backward facing step. The Reynolds 
number based on the length ahead of the step is 1.2 x 10 6 . 
The computational results are presented in Fig. 11 in the 
form of pressure contours and velocity vectors. The pres- 
sure contour plot (Fig. 11a) shows the leading edge shock 
wave created by the boundary layer growth, the expan- 
sion wave at the corner of the step, and the reattachment 
shock. The main and secondary recirculation regions be- 
hind the step are illustrated in the velocity vector plot 
shown in Fig. lib. 

Experimental data consisted of the wall pressure dis- 
tribution and profiles of static pressure at three locations 
downstream of the step. Figure 12a presents the com- 
parison of the wall pressure distribution behind the step. 
The pressure is normalized by the inflow static pressure. 
Figure 12b shows comparisons of the normalized static 
pressure profile at three locations downstream of the step: 
x/h — 0.0537,2.137, and 4.279. These comparisons indi- 
cate good agreement between the CFD and experimental 
results, which verify that the multi-block code can predict 
accurately laminar supersonic flowfields over a step. 

Following this benchmark test case, the code was ap- 
plied to a ram accelerator configuration including a back- 
ward facing step, the schematic of which is shown in 
Fig. 13. A 91 x 45 and 83 x 77 grid is used for the two 
blocks, respectively. The height of the step is 11% of the 
maximum projectile radius, the Mach number of the in- 
cident flow is 4.8, the wall temperature is taken as 600°K 
and the flow is assumed to be laminar. Figure 14 shows 
the nondimensional temperature contours and velocity 
vector field behind the step for the nonreacting solution. 
The reflected shock wave from the tube wall strongly af- 
fects the flow in the recirculation area. Compared with 
the results shown in Fig. 11, it is observed that the flow 
now circulates in the opposite direction. The shock wave 
creates a strong pressure gradient and modifies the cir- 
culation pattern. The reacting flow solution is shown in 
Fig. 15. Reactions start again in the boundary layer be- 
hind the shock impingement point. The combustion is 
then seen to propagate into the recirculation zone. The 
combustion process expands the recirculation zone and 
generates a secondary shock, which is strong enough to 
ignite the main flow above the recirculating zone and a 
detonation wave is established. The detonation wave then 
continues to propagate upstream. Increasing the height of 
the step could perhaps produce better results, but it will 
also considerably diminish the performance of the ram 
accelerator. Therefore, alternative configurations should 
be investigated as a continuation of this work. 


Conclusions 

A CFD code developed for solving the fully cou- 
pled two-dim ensional/axisymmetric Reynolds-averaged 
Navier-Stokes equations including finite rate chemistry 
has been used for studying the flow, combustion and per- 
formance characteristics of the superdetonative ram ac- 
celerator concept. Vector extrapolation methods used in 
combination with the basic iterative scheme show poten- 
tial for significant savings in the overall computational 
work. A diagonalized formulation was tested for the com- 
bustion of a premixed H 2 - air supersonic flow over a com- 
pression corner, and was found to be less efficient than 
the original full Jacobian formulation. 

The results obtained for the ram accelerator configu- 
ration indicate that viscous effects are of primary im- 
portance. Combustion always started in the boundary 
layer, and two types of combustion processes, one stable 
and the other unstable, were established depending on 
the inflow conditions. In the unstable case, a detonation 
wave is formed, which propagates upstream and unstarts 
the ram accelerator. In the stable case, a solution that 
converges to steady-state is obtained, in which the com- 
bustion wave remains stationary with respect to the ram 
accelerator projectile. The unstable combustion mecha- 
nism appears to be caused by a sequence of boundary 
layer separation, upstream propagation of the combus- 
tion process through the boundary layer, and generation 
of secondary shock/detonation waves. The process ap- 
pears to be enhanced by turbulence and by increasing the 
fill pressure. The performance characteristics of a stable 
ram accelerator configuration were computed at various 
Mach numbers. The results indicate that a positive net 
thrust can be obtained even when combustion occurs in 
the nose boundary layer. The possibility of stabilizing the 
detonation wave by means of a backward facing step was 
also investigated. This approach appears to be imprac- 
tical for the configuration studied. A more systematic 
investigation of shock wave/boundary layer interactions 
in premixed combustible gases, and of alternative ram ac- 
celerator configurations, should be carried out as a con- 
tinuation of this work. 
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Table 1: Performance characteristics of the ram accelerator 



Mach number 
M 

Pressure thrust 
F 

Skin friction drag 
D 

Net thrust 

Fnet 

Thrust pressure ratio 
<t>t 

1 

6.7 

3.60 

0.55 

3.05 

17.68 

2 

7.5 

2.43 

0.53 

1.90 

12.41 

3 

8.0 

2.09 

0.45 

1.64 

8.21 


f— Oblique-detonation or 
l shock-deflagration wave 



Figure. 1 Schematic of superdetonative ram accelerator 
mode. 




a) Pressure contours 


(b) Convergence history of Li density residual 


Figure. 2 Supersonic reacting flow past a compression 
corner, (grid: 80 x 50). 
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Figure. 3 Nondimensional temperature T/T a 0 (top 
half), and pressure p/p<x> (bottom half) contours 
showing converged solution for reacting turbulent 
flow (for clarity, the vertical direction is magnified 
by a factor of 2). Mixture: stoichiometric .Hz-air. 
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Figure. 4 Nondimensional pressure distribution along 
the projectile surface and tube wall for various Mach 
numbers. 
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Figure. 8 Mach number as a function of vertical dis- 
tance from the projectile surface, at the x/L = 0.35 
station. 





x/L 


Figure. 5 Variation of y + (at points nearest to the wall), 
skin friction coefficient, and heat transfer coefficient 
along the projectile surface. 
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Figure. 7 Nondimensional temperature T/Too (top 
half), and Mach number (bottom half) contours 
showing solution at several iteration stages. Mixture: 
2H-i + f ?2 + 5Ar, (turbulent flow). 




Page intentionally left blank 



T emperature 

1MTL3C0 3CSTTLD OOCM’TcrCDaCsl "TtDODdN^TLDCO SCMTCO COOCSITT CCCOOCVJ TTUDX1C NTO 

»-’-’-’-rMC\^jcscsjnnr3c , :n r rTT , TrrLrLninLnLn^cDC£)cocgr^rN-r > -rwf^-cDCDcoo335g:CTcnOT 



acvjTrocoafvj^Ttflcoacsj'Tircoaf'j'TCDaoccsT'tflcoaN^cflcoacv 
□□□□□—- - -NN rsicscsir: neon n'^T’T'^r ttlojut in loo ir 


Mach Number 


(b) Detail of shock wave/boundary layer interaction 


(a) Flowfield over a ram accelerator 


Figure. 8 Nondimensional temperature T/T 0 0 , and 
Mach number contours showing converged solution 
for reacting laminar flow. Mixture: stoichiometric 
. 02 -air, (pa, = 1 atm). 
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Figure. 9 Nondimensional temperature T/Tqq, and 
Mach number contours showing laminar flow solution 
at several iteration stages. Mixture: stoichiometric 
. 02 -air, (pa, = 10 atm). 
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Figure. 10 Nondimensional temperature T/T a 0 , and 
Mach number contours showing laminar flow solution 
at several iteration stages. Mixture: 2H2 + C >2 + 5Ar, 
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Figure 11 Pressure contours and velocity vector plot for 
the laminar flow over a backward facing step. 
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(a) Normalized wall pressure distribution 


(b) Normalized pressure profiles 


Figure 12 Comparison of CFD and experimental pres- 
sure distribution for laminar flow over a backward 
facing step. 



Figure 13 Computational grid. 
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a) TOO iter. 



(b) Velocity vector plot behind backward facing step Figure lg Nondimensional tem p era ture T/T^ contours 

showing reacting flow solution at several iteration 
stages. Mixture: 2 H 2 + O 2 + 5^4r, ( p <*, = 1 atm). 


Figure 14 Nonreacting converged solution for laminar 
flow over a ram accelerator. Mixture: 2 H 2 + O 2 + 
5 Ar, (poo = 1 atm). 
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